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(Dated: January 26, 2013) 

We investigate the neutron scattering spectrum in iron pnictides based on the random-phase 
approximation in the five-orbital model with a reahstic superconducting (SC) gap, A — 5meV. 
In the normal state, the neutron spectrum is suppressed by large inelastic quasi-particle (QP) 
scattering rate 7* ~ A. In the fully-gapped s-wave state without sign reversal (3++), a hump-shaped 
enhancement appears in the neutron spectrum just above 2A, since the inelastic QP scattering is 
prohibited by the SC gap. That is, the hump structure is produced by the dissipationless QPs for 
QP energy Ek < 3A. The obtained result is more consistent with experimental spectra, compared 
to the results of our previous paper with A — 50meV. On the other hand, both height and weight of 
the resonance peak in the fully-gapped s-wave states with sign reversal (s±) are much larger than 
those observed in experiments. We conclude that experimentally observed broad spectral peak in 
iron pnictides is created by the present "dissipationless mechanism" in the s++-wave state. 

PACS numbers: 74.20.-z, 74.20.Rp, 78.70.Nx 



I. INTRODUCTION 

Since the discovery of superconductivity in iron pnic- 
tides with high transition temperature {Tc)^, substan- 
tial experimental and theoretical works have been per- 
formed to clarify the mechanism of superconductivity. 
The superconducting (SC) gap in many iron pnictides 
is fully-gapped and band-dependent, as shown by the 
penetration depth measurement^ and the angle-resolved 
photoemission spectroscopy (ARPES)'^'"*. The fully- 
gapped state is also supported by the rapid suppres- 
sion in l/Ti (oc T"; n - 4 - 6) below T^^-'^ . On the 
other hand, P-doped Bal228 and LaFePO^'^^ show the 
nodal line behavior (T-linear dependence) in penetration 
depth measurements. In these compounds, Aig sym- 
metry pairing states with accidental nodes are expected 
theoretically. ^^'^'^ 

In iron pnictides, the intra-orbital nesting of the Fermi 
surface (FS) between the hole- and electron-pockets is 
expected to induce the antiferromagnetic (AF) fluctua- 
tions. Taking this fact into account, the fully-gapped 
sign-reversing s-wave state (s±-wave state) mediated by 
the AF fluctuation had been predicted On the other 
hand, we have demonstrated that the orbital fluctua- 
tion mediated fully-gapped s-wave state without sign- 
reversal (s++-wave state) is realized by the inter-orbital 
nesting, by taking the electron-phonon interaction into 
account. ^^'^^ In the latter scenario, the close relation be- 
tween Tc and the crystal structure revealed by Lee^^, 
e.g., Tc becomes the highest when the As4 cluster is reg- 
ular tetrahedron, is automatically explained^^. More- 
over, the latter scenario is consistent with the large 
SC gap on the z^-orbital band in Bal22 systems^^, ob- 
served by bulk-sensitive laser ARPES measurement^®. 
In addition, the orthorhombic structure transition and 
the corresponding shear modulus softening is well ex- 
plained theoretically^^. In newly discovered K2;Fe2Se2 



with Tc 30K, in which only electron-pockets ex- 
ist, orbital-fluctuation-mediated fully-gapped s++-wave 
state^*^ or spin-fluctuation-mediated nodal d^2_y2-waYe 
state^°'^^ had been predicted theoretically. 

Thus, it is important to clariiy the sign of the SC gap 
via phase-sensitive experiments. Nonmagnetic impurity 
effect offers us useful phase-sensitive information. In iron 
pnictides, the SC state survives against high substitution 
of Fe sites by other element (more than 10%).^^"^^ These 
results support the s+-|_-wave state since the s±-wave 
state is very fragile against impurities, similar to other 
unconventional superconductors. Moreover, impu- 
rity driven crossover from s±-wave state to s+^-wave 
state had been discussed in Refs.^^'^^. 

Another promising method is the neutron scattering 
measurement: As discussed by Monthoux and Scalapino 
in Ref.^®, existence of the resonance peak at a nesting 
wavevector Q is a strong evidence for AF fluctuation 
mediated superconductors with sign reversaP®"'^^. The 
resonance occurs under the condition Wj-os < 2A, where 
cjj-cs is the resonance energy and A is magnitude of the SC 
gap. The sharp and large resonance peak has been ob- 
served in many AF fluctuation mediated unconventional 
superconductors, hke high- Tc cuprates'^^"'^^, CeCoIns'^^, 
and UPd2Al3'^^. The measurements of phonon spectral 
function for \uj\ < 2 A would also be useful. "^^ 

Neutron scattering measurements for iron pnic- 
tides had been performed'^®"*'^ after the theoretical 
predictions^*'^^. Although clear peak structures were 
observed in FeSeo.4Teo.6^^, BaFes-o^Co^.Asz''"'^^ and Ca- 
Fe-Pt-As^3, these weights are much smaller than those in 
high-Tc cuprates and CeCoIns. Moreover, the resonance 
condition cjj-cs < 2A is not surely conflrmed since it is 
difflcult to determine the value of A accurately. For ex- 
ample in BaFei.85Coo.i5As2, ijJrcs is observed as lOmeV 
in neutron scattering measurement.^'^ In this material, 
2 A = Ah + Ac, where Ah(c) denotes magnitude of gap 
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on the hole (electron) pocket. It was estimated as llmeV 
according to an earlier ARPES measurement.^® However, 
Ah+Ae was estimated as 7mcV by a recent measurement 
of the specific heat^''. We also obtain Ah + Ae = 7meV 
from a recent penetration depth measurement in Ref.'**, 
by the linear interpolation for x = 0.14 and x = 0.17. 

In our previous paper^^, we revealed that for A = 
50meV a prominent hump structure free from the reso- 
nance mechanism appears in neutron scattering spectrum 
just above 2A in the S-|-+-wave state. The hump struc- 
ture originates from the dissipationless quasi-particles 
(QPs) free from the inelastic scattering in the SC state. 
Although the broad spectral peak observed in iron pnic- 
tides was naturally reproduced based on the s+-|--wave 
state, rather than the s±-wave state, used model param- 
eters wore not realistic. 

In this paper, we investigate the dynamical spin sus- 
ceptibility X*'('^) Q) based on the five-orbital model^^ for 
both S++- and s±-wave states, by improving the method 
of numerical calculation. Using a realistic parameter 
A = 5meV, the obtained results are more realistic than 
our previous results for A = 50meV.^^ In the normal 
state, x^i'^j Q) is strongly suppressed by the inelastic QP 
damping 7* , which is large due to the strong correlation. 
However, this suppression is released in the SC state since 
the inelastic damping 7* disappears for |a;| < 3A. This 
"dissipationless mechanism" induces a hump-shaped en- 
hancement in x''('*^: Q) in the s++-wave state, just above 
2A till ~ 3A. In the s±-wave state, very high and sharp 
resonance peak appears at cjros < 2A even in the case 
of A = 5mcV. We demonstrate that the broad spectral 
peak observed in iron pnictides is naturally reproduced 
based on the s++-wave state, rather than the s±-wave 
state. 

In Sec. III.C, we comment that Nagai et aL^" fail 
to reproduce the spectral gap in the two-particle Green 
function 2A and are therefore unreliable. In appendix, 
we introduce the similar hump structure of the neutron 
scattering spectrum in CeNiSn much below the Kondo 
temperature Tk- This compound is called Kondo semi- 
conductor since the hybridization gap A opens much be- 
low Tk, while it is an incoherent metal with large in- 
elastic scattering above Tk- This is another example of 
the hump structure by the "dissipationless mechanism" , 
since the inelastic scattering is prohibited by the singlet 
gap A. 

We note that numerical results are improved from re- 
sults in the first version of preprint^^, in which the value 
of 7* (e) for 3 A < e < 4A was incorrect in our previous 
numerical calculation. 



II. FORMULATION 

A. Method of calculation 

Now, we study the 10 x 10 Nambu BCS Hamiltonian 
Hk composed of the five-orbital tight-binding model and 



the band-diagonal SC gap introduced in ref.^^. Then, the 
10 X 10 Green function is given by 



g{iujn, k) 



G{iujn,k) F{iujn,k) 
F'^{iujn,k) -G{-iuJn,k) 
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where a;„ = 7rT(2n 4- 1) is the fermion Matsubara fre- 
quency, G (F) is the 5x5 normal (anomalous) Green 
function, and is the self-energy in the d-orbital basis. 
In this paper, we assume that the magnitude of the SC 
gap is band-independent; jA^^j = A. 

Here, we have to calculate the spin susceptibility as 
function of real frequency. Numerically, it is rather easy 
to use the Matsubara frequency method and the numer- 
ical analytic continuation (pade approximation).'*^''*^ In 
the present study, however, we perform the analytical 
continuation before numerical calculation in order to ob- 
tain more reliable results. The irreducible spin suscepti- 
bility in the singlet SC state is given by^^ 
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[tanh— G|^,3(a;+,fc+)pG^(a;,fe) 
+ tanh ^ {x+,k+)Gt,i, {x, k) 

+ tanh ^FJli^{x+,k+)pi^i^{x, k) 
+ tanh ^pIi^ {x+, k+)F}X [x, 

where a:+ = x-\-uj^ fe+ = fe+q, h — \, - ■ ■ ,5 represents the 
d-orbital, and A (R) represents the advanced (retarded) 
Green function, p^, (a;, k) = {G^, (x, fc) — Gj), (a;, fc))/27ri 

and Pil;^\x, k) = {fIP^{x, k) - f[P^{x, k))/2m are one 
particle spectral functions. Here, we divide y^^'^^^ into 
the "Hermite part" x*^' and "non-Hermite part" yp" , 
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Then, x°' and x"" are expressed as 
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[tanh ^6,^(3 {x+, k+)pf^i^{x, k) 
+ tanh |±pG 3 {x+, k+)Qli^ {x, k) 
+ tanh {x+, k+)p^i^i^ {x, k) 



+ tanh |±pF;^(x+, fc+)ef3|^(x, k)\^) 
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X 



0// 



X [pi%(a;+,fc+)Pu;,(a;,fc) 



tanh- 



where we denote 6^, {x, k) = {Gfl, {x, k) + G^, {x, k))/2-K 



and Q^^^^\x, k) ^ (F^p^ix, k) + F^^}''{x, k))/2n. 



We explain that the non-Hermite part satisfies the re- 
lation x°"(w,9) = for < 2A at T = 0. Now, we 
assume w > 0. pfii^{x,k) = is given for |a;| < A 
since the SC gap opens. Then, in order to satisfy 
both pfi!^ {x,k) 7^ and p^;^(x+,fe+) ^ 0, inequali- 
ties \x^\ > A and > A are required. Moreover, 
[tanh 1^ — tanh ^] ^ only when x+ ■ x < 0. In order 
to obtain the finite value of x^"{u), q) in eq. (5), all three 
inequalities must be satisfied for some x. Considering the 
third inequality, the first two inequalities are restricted 
to 



x+ > A, 
X < -A. 



(6) 
(7) 



They are satisfied for some x only when u > 2A. There- 
fore, x°"{i^,Q) for |a;| > 2A, while X°"('^,9) = for 
|w| < 2A. 

In the present numerical study, we calculate exactly 
X°" using eq. (5), and calculate approximately x"' using 
the Hermite part of eq. (6) in Ref.^^. Using this method, 
we can calculate accurately the imaginary part of the 
spin susceptibility as we will discuss later. 

Then, the spin susceptibility x.'^i^^Q) is given by the 
multiorbital random-phase-approximation (RPA) with 
the intraorbital Coulomb U, interorbital Coulomb U' , 
Hund coupling J, and pair-hopping J'^^: 



1,3 



(8) 



where vertex of spin channel Sf^i^ = U, U', J and 
J' for li = I2 = h = I4, li = ^ I2 = I4 , h = h 
I3 = I4 and h = I4 ^ h = hi respectively. Hereafter, 
we put J = J' = 0.15eV, U' = U - 2 J, and fix the 
electron number as 6.1 (10% electron-doped case). In the 
present model, x''(0,q) takes the maximum value when 
q is the nesting vector Q = (tt, tt/S). Due to the nesting, 
X'(0,Q)/X°(0,Q) « 1/(1 - ast) is enhanced; ast (< 1) 
is the maximum eigenvalue of S^x^^{0, Q) that is called 
the Stoner factor. 

In the following, we prove that the non-Hermite 
part of spin susceptibility Imx''(w,q) = [x'*^(w,q') — 
X''^{uJ,q)\/2i is zero for |w| < 2A at T = 0, 
except at the resonance energy Wres for the s±- 
wave state: The spin susceptibility is expressed as 



OR(A)[^ - S'Ox°^(^^]~^- As explained, x°" = is sat- 



X 

isfied for u < 2A. Then, we obtain x*^ = x"^ = 
X°'[l - S°x°']~^- As a result x*^ - x^^ = for w < 2A 

except when det[l — 5'*'x°'] = 0, which is satisfied at 
uj = cji-cs in the ,s±-wavc state. Thus, if we perform 
the numerical calculation of Eqs. (3)-(8) accurately, 
Imx**(w, Q) = should be satisfied for u) < 2 A. 



B. Inelastic QP damping rate 7* 

In strongly correlated systems, x^i^^l) is renormal- 
ized by the self-energy correction. We phenomeno- 
logically introduce a band-diagonal self-energy as z ■ 
ImS^(e) = i7*(e)l, where z = m/m* is the renormal- 
ization factor. First, we estimate the QP damping in the 
normal state from the experimentally observed conduc- 
tivity. From the Nakano-Kubo formula, the conductiv- 
ity is given by tr = A^^(0)i;^/27(0), where 7(0) = 
7* (0) /z is the "unrenormalized" damping at zero energy, 
and A^i/(0) and are the density of states (DOS) and the 
Fermi velocity of the v-th FS, respectively. Using the five- 
orbital model, we obtain p ps (2.07(0) [me V]) fiQcm for 
the inter-layer spacing c = 6A and p « (2.87(0) [me V]) 
/ificm for c = SA.^^'^'^ In table I, we show the T- 
depcndence of p estimated by fitting the experimental 
data below ~ lOOK^-^"^**, and the inelastic damping 7(0) 
is derived from the theoretical relation between p and 
7. For example, in BaFei.85Coo.i5As2 (c w 6A), the un- 
renormalized inelastic damping 7(0) is estimated as 3.7T, 
which is comparable to that in over-doped cuprates. 

Then, we derive the (e, T)-dependences of the "renor- 
malized" inelastic scattering: In the presence of the 
strong spin and orbital fluctuations, the damping follows 
the approximate relation 7*(e) « b{T + \e\/n) according 
to spin (orbital) fluctuation theories. ^^'^^ According to 
Table I, we obtain b ^ 1.9 in BaFei.85Coo.i5As2 if we 
assume z ~ 0.5. In the present study, we use a larger 
value b = 2.5. Note that the result is not so sensitive to 
the value of b. 

In the present numerical study, we assume more simple 
e-dependence of 7*(e) to simpHfy the analysis, justified in 
calculating Imx* for < |w| < 4A. In the normal state, 
we put 



7*(e) = 7o- 



(9) 



In the SC state at T < T^, 7*(e) = for |e| < 3A (= 
a particle-hole excitation gap (2A) plus a single-particle 
excitation gap (A)), while its functional form approaches 
to that of the normal state for |e| > 3 A. Taking these 
facts into account, we put 



7*(e) = «(e)7: 



(10) 



ll.mrm 



where x*^^^^ 



where (i) a(e) < 1 for |e| < 3A, (ii) a(e) = 1 for |e| > 4A, 
and (iii) linear extrapolation for 3 A < |e| < 4 A; see 
Fig. 1. We have confirmed that the obtained results 
are insensitive to the boundary of |e| (4 A in the present 
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Bai-,K,Fe2As2 [Tc = ilKf^ 


BaFei.s5Coo.i5As2 [Tc = 25K]'^ 


LaFeAsO0.89F0.11 [Tc = 


p{T) - p{Q)[iiQ,cm\ 


~ 23T 


~ 7.3r 


~ 4.6r^ 


7(0) [me V] at T 


~ 12T 


~ 3.7T 




7(0) [me V] at Tc 


~ 37 


~ 7.9 


~ 9.3 



TABLE L p(r) — p(0) and "unrenormalized" inelastic damping at zero energy 7(0) (= 7*(0)/z) estimated by fitting the 
experimental data below ~ lOOK.''^"''''. The unit of T is [meV]. 



case) between (ii) and (iii). Since 7*(e) is an increase 
function of T, 7* at T < Tc should be smaller than 7q. 
Here, we derive the values of 7* and 7g from the relations 
7*(e) - 2. 5(T+|e|/7r), by putting e = 3A = 15meV since 
we are interested in the hump structure around oj ~ 3 A. 
Therefore, we put 7o = 7* (3 A) = 20meV at T = 3meV 
in the normal state. Similarly, we put 7* — 7*(3A) = 
lOmeV at T = in the superconducting state. In the 
s_|__|_(±)-wave state, we put A = 5meV for the two hole- 
pockets and A — (— )5meV for electron-pockets. In the 
numerical calculation, we use 3072 x 3072 fc-meshes and 
a(e)7; = O.SmeV (= O.lA) for |e| < 3A. 



C. Hump structure in Imx* due to dissipationless 

QPs {Ek < 3A) 



(a) particle-hole pair: Ep.\i=Ek'+q+Ek'> 2A 




A 2A 3A 



CO 

FIG. 1: (Color online) (a) Schematic inelastic scattering pro- 
cess in the SC state, by creating a particle-hole excitation 
2A. The realization condition is Ek > 3A. (b) Energy de- 
pendences of the DOS, 7*(aj), and Imx^(i^) in the s++-wave 
state. 



Here, we explain an intuitive reason why the QP is 
"dissipationless" for \uj\ < 3A at zero temperatures. In 
Fig. 1 (a), we show an inelastic scattering process, in 
which a QP at k is scattered to k — q, with exciting a 
particle-hole (p-h) pair {k' + q, k'). Since a QP in the SC 
state cannot exist in the thin shell < A, the particle- 
hole excitation energy i?p_h is always larger than 2A. 
Since the energy of the final state E^^q is also larger 
than A, the inelastic scattering is prohibited when Ef. < 
3A. Thus, the relationship j*{uj) = for \uj\ < 3A is 
obtained. Form this relation, the peak of the DOS at 
a; = A for the isotropic SC gap remains to be sharp. 

Then, the dissipationless QPs in the SC state produce 
the hump-shaped enhancement in the spin spectrum. In 
the normal state, Imx" has no gap structure, and it is 
suppressed by the inelastic QP damping 7* induced by 
the strong correlation. In the SC state as illustrated in 
Fig. 1 (b), Imx*^ has the p-h excitation gap 2 A. Since 
the QP is dissipationless for < 3 A in the SC state, 
the suppression in Imx''(w) is released just above the ex- 
citation gap uj > 2A so as to form a hump structure. 
For this reason, a prominent hump structure appears in 
lmx'{uj, Q) just above 2A till ~ 3A in strongly correlated 
,s++-wave superconductors. 



III. NUMERICAL RESULT 

A. Spin susceptibility at the nesting vector 

Q = (^,7r/8) 

Figure 2 shows the obtained Imx^(a;,Q) at the nest- 
ing vector between the hole- and electron-pockets Q = 
(7r,7r/8): We fix T = ImeV hereafter, since the obtained 
results are insensitive to the temperature for T < 3meV. 
In the normal state with 7q = 20(15)meV, the Stoner fac- 
tor is ast = 0.950(0.959) for U = 1.32eV. In the SC states 
with 7; = lOmeV, ast = 0.956 (0.982) in the s++-wave 
(s^-wave) state for U = 1.32eV. In the s±-wave state, 
ast increases due to the coherence factor. Inversely, ast 
in the s++-wave state decreases due to absence of coher- 
ence factor. As shown in Fig. 2, in the normal state 
with 7g = 20meV, the peak position of Imx^ is about 
20 — 25meV, which is consistent with experimental result 
in BaFei.g5Coo.i5As2'"'. Thus, the value of Imx* in the 
normal state with 7q = 15meV is overestimated. 

A broad hump structure appears in the s+^-wave state 
at w > 2A even in the case of A = 5meV, and its 
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FIG. 2: (Color online) (a) oj-dependence of Imx''(i^,Q) at 
Q = (7r,7r/8) for U — 1.32eV in the s++-wave state (A = 
5meV), as well as in the normal state with 70 = 15, 20meV. 
The hump structure is enhanced by considering the high- 
energy dependence of the SC gap, by introducing the cut- 
off energy AE — 20meV, (b) those in the s±-wave state 
(|A| = 5meV), (c) those for U = 1.33eV in the s++-wave 
state, and (d) those in the s±-wave state. 



this fact into account, we introduce a Gaussian cutoff 
= A''exp{-[el/AE]'^} following Refs.^o-se. We put 
AE = 20meV, which correspond to the Fe ion optical 
phonon frequency ujd 20meV employed in the orbital 
fluctuation theory.^^'^^. 

When the cutoff is applied in the s++-wave state, the 
hump structure becomes more prominent as shown in 
Fig. 2 (a). We confirm that the obvious hump appears 
over the normal state even with 7q — 15meV. The en- 
hancement of hump structure originates from the incre- 
ment of the Stoner factor by introducing the cutoff, from 
ast = 0.956 to 0.965 for U = 1.32eV. 

On the other hand, in the s±-wave state, very high 
and sharp resonance peak appears at Wros < 2A even in 
the case of |A| — 5meV as shown in Fig. 2 (b). This 
result is apparently inconsistent with experimental re- 
sults. In order to explain the experimental result by the 
s-l--wave state, large inhomogeneity would be required, 
although the s±-wave state is fragile against inhomogene- 
ity. The height of the resonance peak exceeds lOOeV"^ 
for 0(0)7* = 0.5meV, while it diverges for a(0) — >■ if 
fc-meshes are fine enough. Imx* is slightly suppressed by 
considering the cut off, AE = 20meV. 

We also study the spectra for both s++- and s±- 
wave states with cutoff for U = 1.33cV: In Fig. 2(c) 
and 2(d), we show the results for the normal state with 
7^ = 15meV (ast = 0.965) and 7^^ = 20meV (ast = 
0.956). We also show results for the s++-wave state 
with 7* = lOmeV (ast = 0.971), and s±-wave state with 
7* = lOmeV (ast = 0.984). 

We note that the effect of multiband on Imx*, which 
was discussed in Ref.^'^, is automatically included in our 
calculation. By increasing U from 1.32eV to 1.33eV, the 
hump structure in the s++-wave state is more enhanced. 
Also, the resonance peak in the s±-wave state develops, 
and uJres shifts to lower energy. 
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overall shape is consistent with experimental results^®' . 
We had neglected the energy-dependence of A in the 
previous study However, in reality, the SC gap A 
will be cut off when the energy of the i^-th band 
measured from the Fermi energy exceeds the character- 
istic energy scale of the pairing interaction. To take 



FIG. 3: (Color online) a;-dependence of Imx" {(jJ , Q) at Q = 
(7r,7r/8) for U = 1.3, 1.32, and 1.33eV in the s++-wave state 
with A = 5meV, and the cutoff energy AE — 20meV, as well 
as normal states with 70 = 20meV for each value of U. 

In Figure 3, we confirm that the hump in the s++- 
wave with A = 5meV state is enhanced as the value of U 
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increases. Thus, the hump becomes prominent as system 
comes close to the AF order. 

In this paper, we have calculated introduced in Eq. 
8. To obtain the value of spin susceptibility x"°" ob- 
served in neutron measurements, we have to take the 
spin magnetic moment (= l/ie) and the factor of spin 
degeneracy. Its z-component is x"°" = 2x''[/J.|eV"^] and 
the transverse spin susceptibility is Xif" = 4x^[Mb6V~^]. 



B. Comparison with our previous method 



had verified that the present "approximated RPA" is re- 
liable in our previous paper In Fig. 1 (b) of Ref.^^, we 
had performed the "exact RPA calculation" for both x°' 
and x°" with A — 400meV, and confirmed that overall 
behavior of Imx'*(cj, Q) is well reproduced by the present 
approximated RPA. 

Here, we comment on the q dependence of Imx'*(w, q) 
around q = Q- In our two-dimensional model^^, it is dif- 
ficult to discuss the q dependence of Imx'^Ci^, q) because 
q dependence of Imx'*(w, q) is drastic even in the normal 
state, which is inconsistent with the neutron scattering 
measurements. 
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FIG. 4: (Color online) Imx"(a;,Q) at Q = (7r,7r/8) for 
U = 1.32eV (a) and U = 1.33eV (b) in the s++-wave state 
(A = 5meV). We show the comparison between the present 
improved method and the previous method in Ref.*^, with 
the cutoff {AE = 20meV). 

In the s++-wave state, Imx^(w,(5) = for |a;| < 2A 
at r = as we discussed in Sec. II. This relation is 
correctly satisfied in the present method if we put a(0) — > 
in Eq. (10). In the present method, we perform the 
numerical calculation of x°" using Eq. (5) exactly. In 
fact, in Fig. 4(a) and (b), we verify that the spectral gap 
of Imx*('i^, Q) is well reproduced in the present method 
with a(0)7* = 0.5meV, demonstrating the superiority of 
the present method to the previous method in Ref.^^. 
In the case of s±-wave state, we obtain Imx*(w,(5) oc 
6{uj — LJ-cQs) for |w| < 2 A if the numerical calculation is 
performed accurately. 

In the present paper, we calculate x'^" in eq. (5) ex- 
actly, while x°' is calculated approximately using eqs. (6) 
and (7) in Ref.^^. We consider this is justified since we 



C. Comparison with Nagai et al.' 
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FIG. 5: (Color online) (a) Imx°(i^,q) in the s++-wave state 
(A = 5meV) for U = 1.32eV and 1.33eV obtained in the 
present study. We also plot the data of Nagai et al.^'^ for U = 
1.375eV, by multiplying 0.39. All resuhs are obtained for the 
cutoff energy AE — 20meV. (b) lmx"{oJ, q) at q = (vr, tt) for 
U — 1.3eV in both s++- and s±-wave states with | A| = 5meV 
and 73 = lOmeV, together with the result in the normal state 
for 7o = 20meV. (c) lmx"{u},q) at q = (tTjTt) in both s++- 
and s±-wave states with |A| — 25meV and 7* = 50meV, 
together with the result in the normal state for 70 = 50meV. 

Recently, Nagai et al.^^ had calculated the neutron 
scattering spectrum using the method proposed in Ref.^^, 
and claimed that (i) hump structure in the s-|__(--wave 
state is smeared when 7* ~ lOmeV and A ~ 5meV com- 
pared to the case of A > 25meV, and (ii) resonance peak 
in the s±-wave state becomes very low and broad. More- 
over, they had also claimed that (iii) one can distinguish 

between the s++-wave and the S-| -wave states from the 

spectrum at q = (7r,7r). 

First, we explain that (i) and (ii) are incorrect state- 
ments based on their inaccurate numerical calculation. 
First, their result fails to reproduce the spectral gap of 
Imx" for id < 2 A as shown in Fig. 5 (a). (One can prove 
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rigorously that Imx* = for a; < 2A at T = 0.) Sec- 
ond, the peak position of the result of Nagai et al.^^ is 
about 2A, while it must be higher energy (~ 3A). In the 
s-i--wave state, the resonance peak should be 5 functional 
structure when a(0) in eq. (10) is enough smaller than 
A. Thus, the low and broad resonance peak of Nagai et 
al.^° is far from the exact behavior of the resonance peak. 
In Nagai's results, fine structures in Imx*(w,q) seem to 
be inappropriately smeared in both s++- and s±-wave 
states. 

Next, we comment on the claim (iii). They pointed 

out the spectrum in the s-|--wavc state with A = 5meV 
are different from that in the s++-wave state with A = 
25meV. Here, we show the results of both s++- and s±- 
wave states in Figs. 5, for (b) |A| = 5meV and (c) 
|A| = 25meV. Since both spectra are almost identical, we 
cannot distinguish between the s++- and s±-wave states 
by the spectrum at the wave vector q = (7r,7r) for the 
same A. This result is reasonable because sign of the 
SC gap is preserved through the (tt, tt) shift for both the 
S++- and s±-wave states. Although claim (iii) is based 
on their numerical result in which the hump of the s++- 
wavc state appears only for A > 25mcV, the prominent 
hump appears in the S-|-+-wave state with A = 5meV 
in our improved numerical results as shown in Fig. 2(a) 
and (c). Thus, we conclude it is impossible to distinguish 
between the s++- and s±-wave states with A = 5meV. 



IV. CONCLUSION 

We have studied the dynamical spin susceptibility 
X^(w,Q) in iron-based superconductors for both s++- 
and s±-wave states, by developing more accurate nu- 
merical method and introducing the high-energy depen- 
dence of the SC gap.**^ In the s+^-wave state, the dissi- 
pationless QPs for < 3 A produce a prominent hump- 
shaped enhancement in x**('^) Q) just above 2A till ~ 3A. 
This "dissipationless mechanism" is unrelated to the res- 
onance. The peak energy of the hump will shift to lower 
energy if we consider the band-dependence and/or the 
anisotropy of the SC gap, as we discussed in Ref.^®. 

On the other hand, in the ,s±-wave state, very high and 
sharp resonance peak appears at uires < 2A. In order to 
explain small and broad peaks observed in Refs.^^'^° as 
the resonance peak in the .s±-wave state, sufficient in- 
homogeneity or small SC volume fraction would be re- 
quired. However, the s±-wave state is fragile against in- 
homogeneity. We concluded that the small and broad 
spectral peak observed in iron pnictides is naturally re- 
produced based on the s++-wave state in the absence of 
inhomogeneity, rather than the s-t-wave state. 

In the Comment on the present paper written by Nagai 
and Kuroki on arXiv^^, the authors repeated their claim 
"smallness of the hump in the s+-j--wave state" based on 
the "old method" that was first developed in Ref.^^. In 
Sec. Ill, however, we actually obtained large hump using 
the "new method" , which is mathematically superior to 



the old method. This discrepancy originates from the 
calculation method as well as the numerical accuracy, 
not from the detail of model parameters, as we discussed 
in our Reply on arXiv^*. 
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Appendix A: Hump structure in the neutron 
inelastic scattering for a Kondo semiconductor 
CeNiSn 

In this paper, we have studied the neutron inelastic 
scattering spectrum in iron pnictide superconductors. In 
the s++-wavc SC state, we confirmed that a large hump 
structure appears just above 2A due to the reduction in 
the inelastic QP scattering 7*, which is the most impor- 
tant finding in this paper. 

Then, a natural question is whether such a hump- 
shaped enhancement by "dissipationless mechanism" is 
universal or not. To answer this question, we discuss 
a Kondo semiconductor CeNiSn. Figure 6 (a) shows 
the neutron inelastic scattering spectrum in CeNiSn at 
q = (0,7r, 0) at low temperatures^^. The observed large 
and broad hump structure in CeNiSn^^'^*^ is very similar 
to that in iron pnictides. CeNiSn is an incoherent metal 
with large inelastic scattering above the Kondo temper- 
ature Tk ~ 30K, while it becomes a semiconductor with 
c-f hybridization gap in the single-particle spectrum (A) 
much below Tk. 

The effective model for the CeNiSn is described as the 
periodic Anderson model (PAM) at half-filhng.^^"^^ Ne- 
glecting the /-orbital degeneracy, the PAM is given as 

^ = E + £/ E /La- + E flrMflfii 

k,a fe,(T i 

k,cr 

where Cka{cj^„) and fkcr{fl„) are annihilation (creation) 
operators for c- and /-electrons, respectively. V is the 
c-f mixing potential, and U is the Coulomb interaction 
for /-electrons. Here, the bandwidth is 2. Afutou and 
Hirashima studied this model at half-filling using the 
dynamical mean-field theory (DMFT) and the quantum 
Monte Carlo (QMC)^^. Hereafter, we introduce their 
numerical results and discuss the energy-dependence of 
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Imx^{i^)- Readers can find more detailed explanations in 
the original paper^^. 



iio 

I 

'S5 

s 



(a) 


1 T 


— 1 — ■ — 1 — 1 — 1 — 1 — 1 — 

= 2K (h k 0) lOIIB ; 


T=10K o[ 


\ ° 
k • \ 












T-^. 







FIG. 6: (a) w-dependence of S{Q,u:) = Imx7(l - e"'^/'^) in 
CeNiSn at various temperatures^^, (b) Single-particle spec- 
trum A{uj) at T = 1/16 for [/ = (solid curve), l(open 
circles), 2(open squares), 3(opon diamonds), 4(solid squares) 
and 5 (solid circles)®^. The inset shows the low frequency part, 
(c) Imx''(w) for f7 = 4 at different temperatures®^. 



Figure 6 (b) shows the obtained single-particle spec- 
trum A{uj). For [/ = 0, the hybridization gap in A(uj) 
is A = 0.62. For J7 = 4, A is rcnormalizcd to 0.35 at 
T=1/16(<Tk), while the gap is smeared out by ther- 
mal fluctuations above Tk^^. At T = 0, inelastic QP 
scattering is suppressed by the hybridization gap, such 
that 7*(w) = for < SA^^ in analogy to Fig. 1 (a). 



Figure 6 (c) shows Imx''(cj) for U = A. In the metal- 
lic state at T = 1/4 Tk), Imx'^(ci;) shows a gap- 
less metallic behavior. In the semiconducting state at 
T = 1/16 (^ Tk), in contrast, it shows a spectral gap 
As and the relation Ag w 2A w 0.7 is recognized. At 
the same time, large hump structure emerges around 
Lu ^ 3 A. Because of the absence of spin resonance 
mechanism, its natural explanation is the rcdiiction in 
the inelastic QP scattering (7*(w) = for |w| < 3A), 
as we discussed in Fig. 1 (b). We must stress the 
hump structure in Fig. 6 (c) is obtained exactly in the 
DMFT, by including both the self-energy and vertex cor- 
rections. Therefore, experimental and theoretical studies 
in CeNiSn strongly support the idea of "hump structure 
in the s+^-wave state" given in Fig. 2 (a), that is ob- 
tained by the RPA by introducing the inelastic QP scat- 
tering 7* (w) phenomenologically. 
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